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Abstract 

Q\ ■ We demonstrate that the Evans-Rosenquist formulation of the optical me- 

chanical analogy, so successful in the application to classical problems, also 
describes the motion of massless particles in the Schwarzschild field of general 
relativity. It is possible to obtain the well known equations for light orbit and 
radar echo delay which account for two exclusive tests of Einstein's field equa- 
tions. Some remarks including suggestions for future work are also added. 
I. INTRODUCTION 

The historical optical-mechanical analogy 1 ' 2 has recently been cast into a 

■ familiar form by Evans and Rosenquist. 3-5 This new formulation, based on 
Fermat's principle, provides an interesting approach that can be profitably uti- 
lized in the solution of many classical problems. The approach works either 
way: The well-known ideas and techniques of classical mechanics can be suc- 

^ 1 cessfully applied to the problem of classical optics and vice versa. Such success 

(3J[)! naturally prompts further enquiry as to whether the applicability of the optical- 

mechanical analogy could be extended even beyond the classical regimes. More 
specifically, a curious student might ask: Can the Evans-Rosenquist (ER) for- 

■ mulation be used to describe the phenomenon of light propagation and massive 
0^ | particle motion in general relativity (GR)? In order to decide this question, let 

us note that it addresses two distinct areas of investigation: One has to first ex- 
amine if the a F = ma" optics of ER does at all describe light propagation in 
the GR scenario. If it does, then the first hurdle is overcome and the question 
of applying the optical analogy to massive particle motion becomes meaningful. 
To what extent this analogy would describe planetary motion in GR is a matter 
| of separate investigation. 

In this paper, we shall examine only the first part of the question that 
relates to the motion of massless particles in GR. We choose Schwarzschild field 
of gravity primarily because of its experimental importance. Besides, on many 
occasions, it is used as an ideal example for illustrating fundamental notions of 
GR. 

The exterior Schwarzschild metric, which is a solution of Einstein's GR field 
equations, has been astoundingly successful in describing various gravitational 
phenomena. Famous experimental tests of GR include the bending of light 
rays, perihelion advance of a planet, radar echo delay, and gravitational red 
shift in the environment of the Schwarzschild gravity field generated by a static, 
spherically symmetric material source like the sun. 6 

In the present approach, the only ingredient of GR to be used is the isotropic 
form of the metric relevant to the gravity field (here the Schwarzschild field). 
No sophisticated mathematics involving four vectors, tensors, Christoffel sym- 
bols, or geodesic equations will be necessary. In addition to being simpler and 
hence understandable by a wide range of readers, the ensuing exposition offers 
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a different window through which the GR events can be visualized. From the 
instructional point of view, it is always useful to look at familiar things from as 
many different perspectives as possible. 

It would be worthwhile to display the relevant formulas of the optical- 
mechanical analogy that we are going to use in our paper. This is done in 
Sec. II. Thereafter, as a necessary step, the Schwarzschild gravity field is por- 
trayed as a refractive optical medium in Sec. III. GR equations for the light orbit 
and the radar echo delay in Schwarzschild gravity are derived in Sees. IV and 
V respectively. Finally, some relevant remarks including a few suggestions for 
future work are made in Sec. VI. 

II. OPTICAL-MECHANICAL ANALOGY 

We need not go into the detailed formulation of the "F = ma" optics of 
ER. Instead, for our present purpose, we reproduce only a table which is self- 
explanatory 3 
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As can be seen, the role of time t is played by the ER stepping parameter a 
having the dimension of length so that ~r' is not a velocity. It is a dimensionless 
quantity. The transition between t and a can be accomplished by using the 
relation 

da — —^dt, (1) 
n z 

where n denotes the refractive index of the optical medium, cq denotes the 
vacuum speed of light. Evans 4 has shown that the stepping parameter can 
be physically interpreted as "optical action". Just as a mechanical particle 
progresses in time along its trajectory, light progresses in optical action along 
its ray. 

The identification 

U(?) = -y (2) 

giving the "potential" U essentially bestows a mechanical character to photon 
motions; the only restriction being that the motion corresponds to mechanics 
at "zero total energy." We can also imagine the possibility that, at least in the 
classical regime, optical and mechanical motions can take place under the same 
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form of " force" /force law. An excellent example is probably the Luneburg lens 
in optics and the harmonic oscillator in mechanics 3 . In order to use the ER 
formulation in the GR regime, it would be necessary to associate a single scalar 
function, the optical refractive index n, with the gravity field under considera- 
tion. The "force" on the massless particles moving in that gravity field can then 
be explicitly obtained 7 via the ER expression grad^^j-^ . 

From now on, we shall use primes (') only to designate a radial coordinate 
(r') and not differentiation with respect to the stepping parameter a. All the 
differentiations will be displayed in full. 

III. GRAVITY FIELD AS A REFRACTIVE MEDIUM 

One might wonder what relationaship could there possibly be between two 
entities apparently as diverse as a gravity field and a refractive optical medium? 
But, indeed, there is one! It was guessed by Einstein and Eddington 8 , formally 
developed by Plebanski and others 9 and utilized in the investigation of of spe- 
cific problems by many physicists. 9,10 We shall only quote the results in a form 
that is easily intelligible. Plebanski 9 has shown that, in a curved spacetime with 
a metric tensor g a p, Maxwell's electromagnetic equations can be rewritten as if 
they were valid in a flat spacetime in which there is an optical medium with a 
constitutive equation. More specifically, with regard to light propagation, the 
gravity field behaves as a refractive optical medium. For example, the gravita- 
tional field exterior to a static, spherically symmetric mass M is equivalent to 
an isotropic, non-homogeneous optical medium with a refractive index n given 
by 10 

n2 ^=( 1 + S\ 1 -^Y^ rn = GM/cl (3) 

where, as usual, G is the gravitational constant, r is the Euclidean radial coor- 
dinate. 

Advanced students who are likely to have a fair grasp on the two field theo- 
ries, Einstein's and Maxwell's, along with the details of algebraic manipulations 
should have no difficulty in pursuing the analysis of Plebanski and de Felice 9 ' 10 . 
However, there is a simpler alternative derivation of Eq.(3), given below, that 
demands only the knowledge of the form of the Schwarzschild exterior met- 
ric. One may also take Eq.(3) at its face value and proceed by treating this 
n(r) as just a given choice of the refractive index in the ER formulation. We 
shall essentially investigate the consequences of such a choice in the succeeding 
sections. 

Consider the exterior Schwarzschild metric in standard coordinates (r', 6, ip, t) 
dg2 = / _ 2m\ ^ _r i _2m\ 1 df n _ /2( ^ 2 + sin2 ^ 2) 



= B(r')cldt 2 - A{r')dr 12 - r' 2 (d9 2 + sin 2 Odip 2 ) (4) 
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where r' > 2m. Redefine the radial coordinate as r' — r(p 1 (r) = r (l + 
so that, if r = u _1 and r' = u'^ 1 , then 



(u) = 1 



mu\ 
~2~) 



(5) 
(6) 



The metric (4) can then be rewritten in the so called isotropic coordinates 

(r, 0, <p, t) as 



ds 2 = n 2 (r)c 2 dt 2 -cj)- 2 (r)[dr 2 +r 2 (de 2 +sm 2 9d<p 2 )} = Q 2 (r)c 2 dt' 2 -(p- 2 {r) |dr>| 
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(7) 

where fi(r) = (l + f) 1 (l - § ) . 

The above form of the metric has a conformally Euclidean spatial part. 11 
Therefore, the isotropic coordinate speed of light c(r) at any arbitrary point in 
the gravitational field, obtained by putting ds 2 = 0, is 



c(r) 



d r 
~dt 



co0(r)ft(r). 



Defining the refractive index as n(r) = co/c(r), we have 

which is precisely the same as Eq.(3) above. It is easy to see that 

4> 2 (u') = 2- 4 [l + (l-2mu') 1/2 l 4 , 
n 2 (u') = 1 - 2mu' 



and 



du' = Q(u)<t>(u)du. 



(8) 
(9) 

(10) 
(11) 
(12) 



All the above relationships will be used throughout the paper. It should be 
mentioned that Sjodin 12 has also derived the expression for 0, Q and n by a 
different method based on Rindler's operational definitions of length and time 
in a gravitational field. 

The refractive index n(r) above corresponds to a radial "force law" having 
a magnitude 

urm dU m(m - 4r)(m + 2r) 5 2m 15m 2 27m 3 _ A . 
"P" — — i Li L_ — I 1 l Ofm ^ 

dr 8(m-2r) 3 r 5 r 2 2r 3 2r 4 1 7 ' 



IV. LIGHT ORBIT EQUATION 
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Within the ER framework, the optical analog of the mechanical zero total 
energy or in short "zero total energy" (see the entry in the table in Sec. II) is 
given by 



d r 
da 



n 2 



= 0. (13) 



We straightaway claim that this very equation is the GR light orbit equation 
in Schwarzschild gravity provided n is given by Eq.(9). In order to see that, 
we proceed as follows. Since there is spherical symmetry in the problem, we 
can, without loss of generality, fix a plane in which the orbiting has to take 
place. It is customary 13 to choose 9 = n/2. Writing out Eq.(13) in plane polar 
coordinates, we have 

^ 2 + ^f^ 2 -» 2 M = o. («) 



da J \da 

Noting that Eq.(13) is the first integral of the "equation of motion" ~~r" ^grad^ 2 ^ . 
Since n does not depend on tp, the ip component of the equation of motion yields 
a conserved quantity called by ER as the "angular momentum" h given by 

h = r 2 ^ = const. (15) 
da 

Eliminating the stepping parameter a, and writing r = u _1 , we have 

2 / du\ 2,-2 



+ -°- < 16 > 

Evans and Roscnquist have solved a number of problems in classical optics/mechanics 
using different choices for n(r). Below we shall discuss, in the context of GR, 
some aspects of the Eqs.(13)-(16): 

(i) As ER have already discussed, the optical quantity h is completely dif- 
ferent from the corresponding conserved mechanical quantity Hq — r 2 dip/dt 
which is proportional to the areal velocity. In the optical formalism, we have 
h = c 2 n 2 r 2 dip/dt so that ho °c nT 2 and dip/dt oc n~ 2 r~ 2 . Hence, with our form 
of n(r), Eq.(9), neither the areal velocity ho nor the angular velocity dip/dt 
remains constant. Of course, for the same force law, the optical and mechanical 
zero energy orbits must have the same form, since o or t are ultimately elim- 
inated. We see that the GR light orbit equation has the form of well known 
classical optics equation 14 , although this interesting fact is not widely noticed. 
The reason for this is that in most text books on GR, the light orbit equation 
is given in a completely different form. 

(ii) Let us obtain the familiar text book form. Using Eqs. (5), (11) and (12) 
in Eq.(16), we get 

0" 2 (1 - 2muT 1 2 + u' 2 r 2 - n 2 h- 2 = 0. (17) 
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Using Eq.(9) for n, we find 



dip 

Differentiating with respect to <p, we get 

, dV 



,/'•' • ( ^ I - 2mu /3 - A"* = 0. (18) 



3mu' 2 = 0. (19) 



This is precisely light orbit equation in Schwarzschild gravity in standard coor- 
dinates. Therefore, one can say that the familiar Eq.(19) represents in disguise 
just the optical analog of the classical zero total energy mechanics. This con- 
clusion will find a further justification in the development of Sec.V. 

(iii) It would be of interest to see how the "angular momentum" h is related 
to the conserved GR quantities E and L, associated with the Killing fields d/dt 
and d/dtp, respectively. This would also provide a relationship between the ER 
stepping parameter a and the geodesic affine parameters used by relativists. 
Integration of GR null geodesic equation gives 



d r 
~~dt 



_» 2 2 r>2 & i. J.-2 -2 C ^' 



= c5ftV, SV—=E, <t>-*r~*-f=L, (20) 
dp dp 



where ds 2 = Xdp 2 , A is a constant and p is a new geodesic affine parameter such 
that ds 2 = =>■ A = 0, dp 2 ^ 0. Eliminating t from Eqs.(20), we get 



Cn 2 E- 2 c 



d~r 2 
dp 



-n 2 = 0. (21) 



If we now define 



dp = c^E-^^da (22) 

then Eq.(13) follows immediately from Eq.(21). We also find that h = L/c E. 
From Eqs.(20 and (22), there also follows the ER relation (1) connecting dt and 
da. Further Eq.(22) implies that ds 2 = Acq 2 E~ 2 4>~ i da 2 , giving the connection 
between a and the affine parameter s. 

(iv) By integrating Eq.(19), one obtains all allowable light orbits 15 .; bound, 
unbound, or even the so-called boomerang orbits 16 . There have been many 
observations confirming the GR predictions of the bending of light rays just 
grazing the sun, the amount being Aip ~ 4m/ Rq, where Rq is the solar radius. 
Interested readers may consult any textbook on GR. We shall only make a 
relevant remark here. M011cr 17 has shown that the bending of light rays is 
due partly to the geometrical curvature of space and partly to the variation of 
light speed in a Newtonian potential. In fact, the ratio of the parts is 50 : 50. 
The GR null trajectory equations can be integrated, once assuming a Euclidean 
space with a variable light speed and again a curved space with a constant light 
speed; both contributing just half the amount of the observed bending. In the 
present approach, on the other hand, we are describing light motion by means 
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of a scalar function n(r). Thus with regard to the light propagation, it looks 
as if spatial curvature and Newtonian potential lost their separate identities 
and merged into an equivalent refractive medium. There is, of course, no point 
in asking which one has a physical reality and which one has not; all these 
are mathematical constructs [like the n(r) here] designed only to interpret our 
physical observations. 

(v) Finally, some words of caution. From the similarity of Eqs.(16) and 
(18), one might be tempted to conclude that in the u' coordinates, n is given by 
n{u') = (1 + 2m/i 2 u' 3 ) 1 / 2 , but that would be incorrect. The correct expression 
for n(u') is obtained by using the fact that n(u) transforms as a scalar. Hence, 
one obtains 18 

n(u') = 4(1 - 2mu')- 1/2 [l + (1 - 2mu') 1 ' 2 ]- 2 . (23) 

Also, it must be understood that, with this n(u'), it is not possible to define 
an isotropic coordinate velocity of light in the vl coordinates. From a direct 
calculation with the metric (4), it will turn out that the coordinate velocities of 
light are different in radial and cross radial directions. 

V. RADAR ECHO DELAY 

Let us now go beyond the geometrical shape of the light orbit and consider 
the dynamics along its path. In other words, we shall derive the GR equations 
of motion involving the time t for the propagation of light rays around a static, 
spherically symmetric gravitating mass M. 

-> 2 2 

Once again we start from the ER "zero energy" equation | i?- — \ = 
and claim that it is the GR equation we are looking for provided n is given by 
Eq.(9). To see this, consider Eqs.(13)-(15) and write 

^-J + h 2 r- 2 -n 2 (r) = 0. (24) 
Utilizing the redefinition u — > u' and the expression for n from Eq.(9), we have 

^ u ,-ifdv!_\ +h 2, 1 _ 2m m '2_ 1 = (25) 
\ da J 

At the position r' [= u^ 1 ] of the closest approach to the gravitating mass M, 
we have dr' /da = => du' /da \ u ' = 0, giving h 2 = (1 — 2mu / )~ 1 u' ~ 2 . Putting 
it in Eq.(25), we find 

^ 4 n- 2 u'- 4 (^Q 2 + u' Q - 2 u' 2 (l - 2mu')- 1 - n- 2 = 0. (26) 

Noting from the metric (4) that B(r') = £! 2 and A(r') = £l~ 2 , and using Eq.(l), 
we finally have 

c-, 2 A(r')B- 2 {r') (^f + r^B" 1 (r' ) B~\r') = 0. (27) 
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This is precisely the textbook form of what is known as the radar echo delay 
equation in Schwarzschild gravity.- and, once again, we see that it is still the 
same ER "zero energy" mechanics, only buried in the (r',t) language. Equation 
(27) is integrated to obtain the time t required by the light signal (in practice, a 
radar signal) to travel from one point of space to another. It is also evident that 
the coordinate speed of light is less than what it would be if the gravitating 
mass were absent. In other words, light is slowed down and the travel time is 
longer. All observers will nonetheless measure a photon's speed to be Cq through 
their positions. In a round trip journey around the mass M, there would be a net 
GR delay in the radar echo reception. This GR prediction has been confirmed 
to a great accuracy by sending radar signals from Earth to Mercury at superior 
conjunction and back. 19 

VI. SOME REMARKS 

The contents of the entire paper vividly demonstrate that light propagation 
in Schwarzschild gravity is indeed describable by the language of "F = ma" 
optics, a shorthand for the optical mechanical analogy. It is remarkable that the 
ER formulation, emerging basically from the investigations of classical problems, 
works equally well also in the GR regime. 

We are concerned with the predictions for light propagation that follow ex- 
clusively from Einstein's field equations. On the other hand, the formula for 
gravitational redshift can be derived from special relativity and the principle 
of equivalence alone. The field equations or their solutions need not be used. 
In that sense, the redshift is not a prediction following exclusively from field 
equations, albeit the effect does follow also from the latter. In the present ap- 
proach, the usual principle of equivalence that equates the effects of gravity and 
artificial forces, is translated into an equivalence of the effects of gravitation 
and refractive medium. Using this idea, the gravitational redshift formula can 
be obtained in its familiar form. We shall present the deductions in a separate 
paper. 

It is evident that the considered approach can deal with a reasonably general 
class of gravity fields where the metric can ce cast into an isotropic form. The 
latter makes it possible to construct a scalar function n which acts as a repre- 
sentative of the gravity field in he matter of light propagation. However, the 
whole range of GR solutions corresponding to different field distributions can 
not be made to correspond to such single scalar functions. At best, a detailed 
constitutive tensor for the equivalent medium may be developed 10 . In this case, 
the challenging task would be to generalize the ER equations appropriately. 
Nevertheless, we can list some immediate future works: 

(i) The motion of a light pulse inside and across the body of a spherical 
star can be tackled quite comfortably. The interior of a spherical star (like 
the sun) is described by the Schwarzschild interior metric. It corresponds to 
a refractive medium with index, say, ri\ while the exterior field corresponds to 
n, Eq.(9). Therefore the whole problem is reduced to one of light propagation 
in a composite media with indices n and n\. At the interface, the matching 
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condition is provided by none other than good old Snell's law. The result will 
provide a theoretical idea about the total path of light rays between the sun's 
core and the Earth. 

(ii) For the sake of completeness, we would expect the U F = ma" optics 
to describe also the motion of massive particles (planetary motion) in the 
Schwarzschild gravity field. To that end, efforts arc underway to extend the 
present treatment. The answer to the second part of the question raised in 
sec. I would then be available. 

Finally, it must be emphasized that there is neither any substitute for nor 
shortcut to to the beauty, generality, and richness of Einstein'e general relativity 
theory. One must eventually grasp all the details of the physics, mathematics, 
and the philosophy of this magnificent never ending edifice. On the other hand, 
the language of the ER optical-mechanical analogy has the power to stimulate 
the interests of a wide cross section of readers who do not have a formal train- 
ing in the sophistications of GR. Those who have the training may, however, 
regard the preceding developments as providing yet another avenue to the same 
experimentally verified tests of light motion in GR. The educational importance 
of such alternative points of view cannot be mistaken. 
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